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RESEARCHES IN THE LUNAR THEORY,* 
By Gr. W. Hill, KyacJc Turnpike, N. Y. 

When we consider how we may best contribute to the advancement of 
this much -treated subject, we cannot fail to notice that the great majority of 
writers on it have had before them, as their ultimate aim, the construction 
of Tables: that is they have viewed the problem from the stand-point of prac- 
tical astronomy rather than of mathematics. It is on this account that we 
find such a restricted choice of variables to express the position of the moon, 
and of parameters, in terms of which to express the coefficients of the peri- 
odic terms. Again, their object compelling them to go over the whole field, 
they have neglected to notice many minor points of great interest to the 
mathematician, simply because the knowledge of them was unnecessary for 
the formation of Tables. But the developments having now been carried 
extremely far, without completely satisfying all desires, one is led to ask 
whether such modifications cannot be made in the processes of integration, 
and such coordinates and parameters adopted, that a much nearer approach 
may be had to the law of the series, and, at the same time, their convergence 
augmented. 

Now, as to, choice of coordinates, it is known that, in the elliptic theory, 
the rectangular coordinates of a planet, relative to the central body, the axes 
being parallel to the axes of the ellipse described, can be developed, in terms 
of the time, in the following series, 

sz=:aX — J \ cos iff, 

i= — CO 1i ~2~ 

i= + CO ]_ («— 1) 

y — l^ _ J ie sin iff, 

i= — CO 1 "2 - 

a and i being the semi-axes of the ellipse, e the eccentricity, g the mean 
anomaly and, for positive values of i, the Besselian function (in Hansen's 
notation) 

:: 1.2. ..i L l.(i 4-1) + 1.2. ft +1) (i 4- 2) —•'"]' 
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6 Hill, Researches in the Lunar Theory. 

while, for negative values of i, we have 

A —A ' 

and, for the special case of i ■=. 0, we put the indeterminate 

% ~ 2 e - 

Here the law of the series is manifest, and the approximation can easily be 
carried as far as we wish. But the longitude and latitude, variables em- 
ployed by nearly all the lunar- theorists, are far from having such simple 
expressions ; in fact, their coefficients cannot be expressed finitely in terms of 
Besselian functions. And if this is true in the elliptic theory, how much 
more likely is a similar thing to be true when the complexity of the problem 
is increased by the consideration of disturbing forces? We are then justified 
in thinking that the coefficients of the periodic terms in the development of 
rectangular or quasi-rectangular coordinates are less complex functions of 
their parameters than those of polar coordinates. There is also another 
advantage in employing coordinates of the former kind ; the differential 
equations are expressed in purely algebraic forms ; while, with the latter, 
circular functions immediately present themselves. For these reasons I have 
not hesitated to substitute rectangular for polar coordinates. 

Again, as to parameters, all those who have given literal developments, 
Laplace setting the example, have used the parameter tn, the ratio of the 
sidereal month to the sidereal year. But a slight examination, even, of the 
results obtained, ought to convince any one that this is a most unfortunate 
selection in regard to convergence. Yet nothing seems to render this para- 
meter at all desirable, indeed, the ratio of the synodic month to the sidereal 
year would appear to be more naturally suggested as a parameter. Some 
instances of slow convergence with the parameter m may be given from 
Delaunay's Lunar Theory: the development of the principal part of the 

coefficient of the evection in longitude begins with the term -j me, and ends 

with the term i 5 28823808o m&e ' a S' a ^ n » ^ n the principal part of the coefficient 
of the inequality whose argument is the difference of the mean anomalies 

21 

of the sun and moon, we find, at the beginning, the term -j mee 1 , and, at the 

end, the term — 8538944 — ^etf. It is probable that, by the adoption of some 
function of m as a parameter in place of this quantity, whose numerical 
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value, in the case of our moon, should not greatly exceed that of m, the 
foregoing large numerical coefficients might be very much diminished. And 
nothing compels us to use the same parameter throughout ; one may be used 
in one class of inequalities, another in another, as may prove most advan- 
tageous. It is known what rapid convergence has been obtained in the series 
giving the values of logarithms, circular and elliptic functions, by simply 
adopting new parameters. Similar transformations, with like effects, are, 
perhaps, possible in the coefficients of the lunar inequalities. However, as 
far as my experience goes, no useful results are obtained by experimenting 
with the present known developments; in every case it seems the proper 
parameter must be deduced from a priori considerations furnished in the 
course of the integration. 

"With regard to the form of the differential equations to be employed, 
although Delaunay's method is very elegant, and, perhaps, as short as any, 
when one wishes to go over the whole ground of the lunar theory, it is subject 
to some disadvantages when the attention is restricted to a certain class of 
lunar inequalites. Thus, do we wish to get only the inequalities whose coeffi- 
cients depend solely on to, we are yet compelled to develop. the disturbing 
function R to all powers of e. Again, the method of integrating by unde- 
termined coefficients is most likely to give us the nearest approach to the 
law of the series ; and, in this method, it is as easy to integrate a differential 
equation of the second order as one of the first, while the labor is increased 
by augmenting the number of variables and equations. But Delaunay's 
method doubles the number of variables in order that the differential 
equations may be all of the first order. Hence, in this disquisition, I have 
preferred to use the equations expressed in terms of the coordinates, rather 
than those in terms of the elements ; and, in general, always to diminish the 
number of unknown quantities and equations by augmenting the order of the 
latter. In this way the labor of making a preliminary development of R in 
terms of the elliptic elements is avoided. 

In the present memoir I propose, dividing the periodic developments of 
the lunar coordinates into classes of terms, after the manner of Euler in his 
last Lunar Theory,* to treat the following five classes of inequalities : — 

1. Those which depend only on the ratio of the mean motions of the 
sun and moon. 



* Theoria Motuum Lunce, nova meikodo pertraciata. Petropoli, 177$. 
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2. Those which are proportional to the lunar eccentricity. 

3. Those which are proportional to the sine of the lunar inclination. 

4. Those which are proportional to the solar eccentricity. 

5. Those which are proportional to the solar parallax. 

A general method will also be given by which these investigations may be 
extended so as to cover the whole ground of the lunar theory. My methods 
have the advantage, which is not possessed by Delaunay's, that they adapt 
themselves equally to a special numerical computation- of the coefficients, 
or to a general literal development. The application of both procedures 
will be given in each of the just mentioned five classes of inequalities, so 
that it will be possible to compare them. 

I regret that, on account of the difficulty of the subject and the length 
of the investigations it seems to require, I have been obliged to pass over 
the important questions of the limits between which the series are con- 
vergent, and .of the determination of superior limits to the errors committed 
in stopping short at definite points. There cannot be a reasonable doubt 
that, in all cases, where we are compelled to employ infinite series in the 
solution of a problem, analysis is capable of being perfected to the point of 
showing us within what limits our solution is legitimate, and also of giving 
ixs a limit which its error cannot surpass. When the coordinates are devel- 
oped in ascending powers of the time, or in ascending powers of a parameter 
attached as a multiplier to the disturbing forces, certain investigations of 
Cauchy afford us the means of replying to these questions. But when, for 
powers of the time, are substituted circular functions of it, and the coefficients 
of these are expanded in powers and products of certain parameters pro- 
duced from the combination of the masses with certain of the arbitrary 
constants introduced by integration, it does not appear that anything in the 
writings of Oauchy will help us to the conditions of convergence. 



Chapter I. 
Differential Equations. — Properties of motion derived from JacoWs integral. 

We set aside the action of the planets and the influence of the figures 
of the sun, earth and moon, together with the action of the last upon the 
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sun, as also the product of the solar disturbing force on the moon by the 
small fraction obtained from dividing the mass of the earth by the mass of 
the sun. These are the same restrictions as those which Delaunay has 
imposed on his Lunar Theory contained in Vols. XXVIII and XXIX of the 
Memoirs of the Paris Academy of Sciences. The motion of the sun, about 
the earth, is then in accordance with the elliptic theory, and the ecliptic is 
a fixed plane. 

Let us take a system of rectangular axes, having its origin at the centre 
of gravity of the earth, the axis of x being constantly directed toward the 
centre of the sun, the axis of y toward a point in the ecliptic 90° in advance 
of the sun, and the axis of z perpendicular to the ecliptic. In addition, we 
adopt the following notation : — 

f = the distance of the sun from the earth ; 

W zr the sun's longitude ; 

(i = the sum of the masses of the earth and moon, measured by the 
velocity these masses produce by their action, in a unit of time, 
and at the unit of distance; 

m! = the mass of the sun, measured in the same way ; 

n' = the mean angular velocity of the sun about the earth ; 

a' = the sun's mean distance from the earth. 
In accordance with one of the above-mentioned restrictions we have the 
equation : — 

The axes of x and y having a velocity of rotation in their plane, equal 

to — -, it is evident that the square of the velocity of the moon, relative to the 
at 

earth's centre, has for expression, in terms of the adopted coordinates, 

* d^'l 2 .fdy dX'-i 2 dz 2 

2 1 - lit ~ y m\ + \jt + x dt] + df 

da? + dtf + dz 2 . dl! xdy — ydx d7J 2 . „ 

= — dV + 2 dt - JL W L ~ + IF V+fr 

The potential function, in terms of the same coordinates, is 

li n a «' 3 n' 2 a' s 

~ V (x 2 + y 2 + z 2 ) + V[(^— W+f + ^\ ^~ X ' 

If the second radical in this expression is expanded in a series proceeding 

• VI ft 

according to descending powers of r', and the first term — - omitted, since it 
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disappears in all differentiations with respect to the moon's coordinates, the 
following expression is obtained 



a = 



V(* 2 + y 2 + z 2 ) 



+ «' 2 ^[}' 2 -*(3' 2 + z 2 )] 

+ |l ~ s \f ~ ^ if + *) + * (f + z 2 Y] 
+ ~ \f - f& (y 2 + * 2 ) + * * (y» + , 2 ) 2 ] 



+ 

Since the differential equations of motion are of the form 

dt ' d. d -* d$~~ dq> ' 

4> denoting, in succession, each of the variables which define the position of 
the moon, it is plain that the term 

may be removed from T and added to £l ; and these modified quantities may 
be denoted by the symbols T' and Of. Then these equations may be written 
thus: 

<Px c> dK dy dft! _ d£V 

lf~" It dt ~ If y ~~ ~dx' 

d 2 y , d%' dx d*X __ dQf 

If + dt Tt + If X ,== ~dy' 

dh dD,' 

If = Hz' 

When we wish to restrict our attention to the lunar inequalities which 
are independent of the solar parallax, all the terms, in the last expression of 
fl, which are divided by f A , r' 5 , r'\ &c, may be omitted. In this case it will 
be seen that all the terms, introduced into the differential equations by the 
solar action, are linear in form, with variable, but known coefficients, since 

— , — --„ and — are known functions of t. 
dt dt- ¥* 
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When all the inequalities, involving the solar eccentricity, are neglected, 
the equations admit an integral in finite terms. For, in this case, we have 

dfr d 2 X 

dt= n > !f = ' f = ^ 

and Cl' does not explicitly contain t; hence, multiplying the equations sever- 
ally by the factors dx, dy, and dz, and adding the products, both members 
of the resulting equation are exact differentials. Integrating this equation, 

we have 

dx 2 + dij 1 -f dz 2 

9Jt2 — " + a constant. 

This integral equation appears to have been first obtained by Jaeobi.* 
As it will be frequently referred to in what follows, I shall take the liberty 
of calling it Jacobi's integral. 

If we take two imaginary variables 



u = x + V (y 2 + z 2 ) V — 1, 
s — x — V (y 2 + z 2 )- V — 1, 



CI has the following simple expression, being a function of two variables only, 

_- « w'V 3 «'V 3 . . 

CI = ~~p= H —J- r -j- r -jTJT {U + S). 



n a a* 



If this is expanded in descending powers of f, and,. as before, the term 
omitted, 

CI = -£= -f n' 2 -J; l~i u 2 + 1 us + I s 2 ] 



[t 8 A « 4 + & ^ + A mV + 3% MS ' + t 8 A « 4 ] 



. a' 2 z 5 

+ 

The additional variable, necessary to complete the definition of the moon's 
position, may be so taken that the expression of T may be simplified as much 
as possible. This expression may be written 

~ m duds a ( ydz — zdy ) 2 , „ d% xdy — ydx , d?f , 2 \ 

— ~W {u—sfdt 2 Tt di "*" df { } ' 

* Comptes Rendus de I'Acudemie des Sciences de Paris. Tom. iii, p. 59. 
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There does not seem to be any function of x, y and z, which, adopted as a 
new variable to accompany u and s, would reduce this to a very simple form. 
However, when we are engaged in determining the inequalities independent 
of the inclination of the lunar orbit, this transformation will be useful to us. 
For, in this case, z = 0, and the values of u and s become 

u = x + y V — 1, 
s = x — y V — 1, 
and T is given by the equation 

2 rp duds d'Af uds — sdu . d% a 

~~df It W~ ~ + W ws ' 

Although n is expressed most simply by the systems of coordinates we 

have just employed, the integration of the differential equations will be 

easier, if we suppose that the axes of x and y have a constant instead of a 

variable velocity of rotation, the axis of x being made to pass through the 

mean position of the sun instead of the true. To obtain the expression for 

d"X' 
T correspondent to this supposition, it is necessary only to write n' for — in 

at 
the former values. As for fl, it can be written thus 

o -- (L _l fl-V 8 __ A" a 

where 

r 2 — x 2 + y 2 -f z 2 = the square of the moon's radius vector ; 

S = x cos v + y sin v ; 

v = the solar equation of , the centre. 
This function being expanded in a series of descending powers of r', as before, 
we have 

DJ = H + 4 n' 2 (x 2 + f) 

+ 



Hill, Researches in the Lunar Theory. 13 

And the corresponding differential equations are 

d 2 x c, ,dy dDf 

dt 2 dt dx 

d 2 y 9 , dx d£V 

W + dt~~dy~' 

<Pz _d& 

If ~~ dz ' 

Thus much in reference to the equations under as general a form as we 
shall have occasion for in the present disquisition. We shall now suppose 
that they are reduced to as restricted a form as is possible without their 
becoming the equations of the elliptic theory ; that is, we shall assume that 
the solar parallax and eccentricity and the lunar inclination vanish. With 
these simplifications, in the first system of coordinates, 

rp, dx 2 -\- dy* . , xdy — ydx 

1 - w' + ,,— #~ ■ 

rp, duds n' uds — sdu 

~ = 2di 2 ' 2 dt ' 

11' = -£=+ !»"(« + a) 2 . 

s/us 

And the differential equations, correspondent, are, in the first case, 

g — 2n' d l + [4 - 3»' 2 ] x = 0, 
dt 2 dt Lr J 



and, in the second, 



d *y A-2rt dx A- £ v — 



and, in the second, 



The Jacobian integral has severally the expressions 

dx -\- dy [i . , 2 2 p 

2dr r 

tt = -4, + f »' 2 (« + s) 2 — C. 
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The terms — 2n'-/, 2n' -— , &c, have been introduced into the equations 
at at 

by making the axes of coordinates movable ; but since the putting of n' = 
makes the solar disturbing force vanish, there is no inconsistency in attribu- 
ting them to the solar action. Then, in the case of the vanishing of this 
action, we have the equations of ordinary elliptic motion 



dt 2 r 



, + 5y=o- 



Thus, in the restricted case we consider, all the terms, added to the differ- 
ential equations of motion by the solar action, are linear in form and have 
constant coefficients. This, and the circumstance that t does not explicitly 
appear in the equations, are two advantages which are due to the particular 

selection of the variables x and y. If -^ were constant, the equations would 

be linear with constant coefficients and easily integrable. 

The constants ^ and n' can be made to disappear from the differential 
equations, if, instead of leaving the units of length and time arbitrary, we 
assume them so that (i = 1, and n' = 1. The new unit of length, will then 

be equal to f^ units of the previous measurement. The equations, thus 
simplified, are 

■*!+G-»].=* 



d 2 x 
It 2 



with their integral 



dx 2 + dtf = 2_ + 3 ^_2(?. 
dt 2 r 



It will be perceived that, in this way, we make the differential equations 
absolutely the same for all cases of the satellite problem. 

O.0C 

Let us put p = — , then 



dt 



fy _ 9 T 2 _l ?ua ___ cyr< »1 i 

dt 



2 [^ + 3x 2 -2C- P f-[l-3]x. 
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Or, by making y the independent variable, 

dx p 

ri _ 3 ~i ^ 

^l — 2 - 



^ r 2 _i_^_^ yi* 



[_ + 3a --2(7-^] i 



The problem is then reduced to the integration of two differential equations 
of the first order. Were this accomplished, and p eliminated from the two 
integral equations, we should have the equation of the orbit. If we put 

W= 2x+[- + Sx 2 — 2C— p 2 ] \ 

the differential equations can be written in the canonical form, 

dx _ dW 

dy dp ' 

dp _ dW 

dy dx 

It may be worth while to notice also the single partial differential 
equation, to the integration of which our problem can be reduced. Returning 
to the arbitrary linear and temporal units, and, for convenience, reversing 
the sign of C, if a function of x and y can be found satisfying the partial 
differential equation 

|"4^ + n'yT + [*I — n'xT = , 2 / ^ + M V + 2C, 
L dx * J L dy J V (x 2 + f) ' 

and involving a single arbitrary constant h, distinct from. that which can be 
joined to it by addition, the intermediate integrals of the problem will be 
dx dV , dy dV 

"dt~^ + n ^ Tt~~dJ~ nx ' 

and the final integrals 

dV dV 4 , 

-dh =a > *? = '+*' 

a and c being two additional arbitrary constants. The truth of this will be 
evident if we differentiate the four integral equations with respect to t and 
compare severally the results with the partial differential coefficients of the 
partial differential equation with respect to x, y, h and C. 
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Although, in this mariner, the problem seems reduced to its briefest 
terms, yet, when we essay to solve it, setting out with this partial differential 
equation, we are led to more complex expressions than would be expected. It 
would be advisable, in this method of proceeding, to substitute polar for 
rectangular coordinates, or to put 

x = r cos <p, y = r sin (p. 

The partial differential equation, thus transformed, is 

4t + ["- d -^— % ' r T = — + § n> v + 2 °+ § ***? c ° s %$■ 

dr \-r d<p A r T 

This would have to be integrated by successive approximations, and it is 
found that this method, which, at first sight, seems likely to afford a briefer 
solution of the problem, because but one unknown quantity was to be deter- 
mined, and this free from the variable t, and involving only half of the number 
of arbitrary constants introduced by integration, when developed, leads to as 
complex operations as the older methods, and, moreover, has the disadvantage 
of giving results which need prolix transformations before the coordinates can 
be exhibited in terms of the time. 

Although we shall make no use of the equations in terms of polar 
coordinates, they may be given here, for the sake of some special properties 
they possess in this form. They are 

r 4S — r2 !§ — 2»'r 2 § + ^— 3»'VcosV = 0, 
at 2 at 2 dt r 

Tt \f ^dt + *) 1 + f n ' V sin 2tp = °' 

with their integral 

' -*- = -£ + 3n' 2 r 2 cos 2 * — 2C, 

dr r 

By dividing the second of the differential equations by r ! , the variables are 
separated, and X denoting the longitude of the moon, we have 

K — in' 2 J 5* — at 

r = — _= s at 

X dt 
K being a constant. Thus, after the longitude is determined in terms of t, 
the radius vector is obtained by a quadrature. But it can also be found, 
without the necessity of an integration, by dividing the' integral by r 2 and 
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1 fir' 1 
then eliminating the term — — , by means of its value derived from the 

r 2 dtf 

second differential equation ; in this way we get 

r 72 

C 



ft 

«,3 



r - 



"|| + i«' 2 sin2<?>" 



It +n 



' d¥ 



+ i^ l2 - — f^cos> 



As we desire to make constant numerical application of the general 
theory, established in what follows, to the particular case of the moon, Ave 
delay here, for a moment, to obtain the numerical values of the three 
constants p, ri and C. The value of fi may be derived either from the 
observed value of the constant of lunar parallax combined with the mean 
angular motion of the moon, or from the intensity of gravity at the earth's 
surface and the ratio Of the moon's mass to that of the earth. We will adopt 
the latter procedure. The value of gravity at the equator, g = 9.779741 
metres, the unit of time being the mean solar second. We propose, however, 
to take the mean solar clay as the unit of time, and the equatorial radius of 
the earth as the linear unit. This number must then be multiplied by 

iJnm\& ' (6377397.15 metres is Bessel's value of the equatorial radius.) 
Moreover, the theory of the earth's figure shows that, in order to obtain the 
attractive force of the earth's mass, considered as concentrated at its centre 
of gravity, a second multiplication must be made by the factor 1.001818356. 
With our units then this force is represented by the number 11468.338 : and 

the moon's mass being taken at 622 _ 7 of the earth's, her attractive force is 
represented by the number 140.676. Consequently 

p = 11609.011. 

The sidereal mean motion of the sun in a Julian year is 1295977".41516, 

whence 

n' = 0.017202124. 

The value of O might be obtained from the observed values of the moon's 
coordinates and their rates of variation at any time. However, as the eccen- 
tricity of the earth's orbit is not zero, C obtained in this manner would be 
found to undergo slight variations. The mean of all the Vcxlues obtained 
in a long series of observations might be adopted as the proper value of this 
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quantity when regarded as constant. But it . is much easier to derive it 
approximately from the series 

2C= fan)* [i + 2m — « m 2 — m 3 — ijg m 4 — ^m 5 - f§*§ m 6 — »» m 7 ], 
which will be established in the following chapter. Here n denotes the 

/ 

moon's sidereal mean motion, and in is put for ■. In this formula the 

n — n 

terms which involve the squares of the lunar eccentricity and inclination and 
of the solar parallax are neglected ; this, however, is not of great moment, 
as, being multipled by at least m 2 , they are of the fourth order with respect 
to smallness. The observations give n = 0.22997085, hence 

C= 111.18883. 

If it is proposed to assume the units of time and length so that [i and 
n may both be unity, it will be found that the first is equal to 58.13236 
mean solar days, and the second to 339.7898 equatorial radii of the earth. 
The corresponding value of G is 3.254440. 

Let lis now notice some of the properties of motion which can be derived 
from Jacobi's integral. This integral gives the square of the velocity relative 
to the moving axes of coordinates ; and, as this square is necessarily positive, 
the putting it equal to zero gives the equation of the surface which separates 
those portions of space, in which the velocity is real, from those in which 
it is imaginary. This equation is, in its most general form, 

which is seen to be of the 16th degree. As y and z enter it only in even 
powers, the surface is symmetrically situated with respect to the planes of xy 
and xz. The left member is necessarily positive, (the folds of the surface, for 
which either or both the radicals receive negative values, are excluded from 
consideration), hence the surface is inclosed within the cylinder whose axis 
passes through the centre of the sun perpendicularly to the ecliptic, and 
whose trace on this plane is a circle of the radius 

As, in general, the second term of the quantity, under the radical sign, is 
much smaller than the first, this radius is, quite approximately V3 a . Thus, 

in the case of our moon, assuming — 7 = sin 8". 848, we have this radius = 
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V3.001383 a . It is evident . that, for all points without this cylinder, the 
velocity is real ; and as, for large values of z, whether positive or negative, 
the left member of the equation becomes very small, it is plain that the 
cylinder is asymptotic to the surface. Every right line, p'erpendicular to the 
ecliptic, intersects the surface not more than twice, at equal distances from 
this plane, once above and once below. 

Let us now find the trace of the surface on the plane of xy. Putting 
p for the distance of a point on this trace from the centre of the sun, 

f = { a'-xf + y\ 
and it is evident that the cubic equation, 

will give the limits between which the values of p can oscillate. If C is 
negative, this equation has but one real root which is negative ; consequently, 
in this case, the surface has no intersection with the plane of xy. But, in all 
the satellite systems we know, C is positive, and . this condition is probably 
necessary to insure stability. Hence we shall restrict our attention to the 
case where C is positive. Then all the roots of the equation are real, and 
two are positive. It is between the latter roots that p must always be found. 
To compute them, we derive the auxiliary angle d from the formula 

S me = [i + ^<L]~\ 

or, since d differs but little from 90°, with more readiness from 

9 ^_ Pi -U a -9— 4- -*- ——~\ 
n' 2 a' 2 L + 3 n' 2 a' 2 "*" * 7 w'V 4 J 



COS 2 



t 1 + f *V J 



C 
' is a small quantity, with sufficient approximation from 



/ 9 / 9 

n z a i 



4 



a 



cos d = 



1 _i_ a ^_ 

The two roots are then 

C \ ■ d 



p 1 = 2«v(l + i^ 2 ) sin 3 

p 2 = 2A/(l + t^>m(60°-|) 



20 
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The trace of the surface on the plane of xy is then wholly comprised in the 
annular space between the two circles described from the centre of the sun as 
centre with the radii p x and p 2 . Moreover, as in most satellite systems we 

haVG A~ 3 GqUal t0 a V6ry Sma11 fraction ' ( for our moon ^v~3 = 3229302)' 
it is plain that, for points whose distance from the earth is comparable with 
their distance from the sun, the trace is approximately coincident with these 

circles. For the term ■!—, in the equation, may then be neglected in com- 

r 

parison with the other terms. 

In the case of our moon there is found 

= 87° 52' 11".53. 
and hence 

p x = 22815.15. p 2 = 23816.09, 

and, if r and p are regarded as the variables defining the position of a point 
in the plane xy, the following table gives some corresponding values of these 
quantities, for each of the two branches of the trace approximating severally 
to the two circles. 



r. 

433.3251 


P- 

22818.69 


r. 

439.1922 


P- 
23151.81 


450 


22816.11 


450 


23153.31 


500 


22869.68 


500 


23160.04 


600 


22860.13 


600 


23169.85 


1000 


22841.59 


1000 


23188.81 


0000 


22811.10 


10000 


23813.43 


.6121.10 


22815.68 


41121.55 


23815.53 



The first and last values correspond to the four points where the curves 
intersect the axis of x on the hither and thither side of the sun. It will be 
seen that the approximation of the branches to the circles is quite close, 
except in the vicinity of the earth, where there is a slight protruding away 
from them. 

In addition to these two branches of the trace, there is, in the case where 
C exceeds a certain limit, a third closed one about the origin much smaller 
than the former. As the coordinates of points in this branch are small 
fractions of a', its equation may bo written, quite approximately, 



(i 



— C—%n'\x\ 
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It intersects the axis of y at a distance from the origin very nearly 

v — * 

and the axis of x at points whose coordinates are the smallest (without regard 
to sign) roots of the equations 



„'2,/ 3 



H + !L^_ = c+ I »'V a — * «' 2 («' — x) 2 , 

_ £ + *?*L = c+ * »' V*— I »' 2 («' — x) 2 

x a — x 

For the moon these quantities have the values 

# = 104.408, ^ = —109.655, x 2 = +109.694. 

This branch then does not differ much from a circle having its centre at the 
origin, more closely it approximates to the ellipse whose major axis = x % — x t , 
and minor axis = 2y . 

The value of the coordinate z, for the single intersection of the surface 
with the axis of z above the plane of xy, is given by the single positive root 
of the equation 

(i n a ~ , ,2 ,2 

For the moon the numerical value of this root is 

z = 102.956. 
The intersection of the surface with the perpendicular to the plane of xy 
passing through the centre of the sun is, in like manner, given by the equation 

+ — = C+ I n'*a>\ 



V(«' 2 + 2 2 ) Z 

having but a single positive root, which is nearly 



la! 



1 + ^ — j-72 

n a 



or, with less exactitude, 

From these investigations it is possible to get a tolerably clear idea of 
the form of this surface. When C exceeds a certain limit, it consists of three 
separate folds. The first being quite small, relatively to the other two, is 
closed, surrounds the earth and somewhat resembles an ellipsoid whose axes 
have been given above. The second is also closed, but surrounds the sun, 
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and has approximately the form of an ellipsoid of revolution, the semiaxis in 
the plane of the ecliptic being somewhat less than «', and the semiaxis of 
revolution perpendicular to the ecliptic and passing through the sun being 
about two-thirds of this. This fold has a protuberance in the portion neigh- 
boring the earth. The third fold is not closed, but is asymptotic to the 
cylinder mentioned at the beginning of the investigation of the surface. Like 
the second, it also is nearly of revolution about an axis passing through the 
centre of the sun and perpendicular to the ecliptic. The radius of its trace 
on the ecliptic is about as much greater than a', as the radius of the trace 
of the second fold falls short of that quantity. The fold has a protuberance 
in the portion neighboring the earth, and which projects towards this body. 
The whole fold resembles a cylinder bent inwards in a zone neighboring the 
ecliptic. 

What modifications take place in these folds when the constants involved 
in the equation of the surface are made to vary, will be clearly seen from 
the following exposition. Let us, for brevity, put 

n' 2 a' 2 
and, for the moment, adopt a', the distance of the earth from the sun, as the 
linear unit, and transfer the origin to the centre of the sun, and moreover put 

a 

*v = — - — . 
r n' 2 a' s 

Then the intersections of the surface, with the axis of x, will be given by the 

two roots of the equation 

a* — X s — hx 2 +(h + 2 — 2y)x — 2 = 0, 
which lie between the limits and 1 ; by the two roots of 

a* — a? — hx 2 + (h + 2 + 2y)x — 2 = 0, 
which lie between 1 and s/h ; and by the two roots of 

x* — x 3 — to?. + (h — 2 — 2y) x + 2 = 0, 
which lie between and — sfh. 

Hence, if C diminishes so much that the first of these three equations 
has the two roots, lying between the mentioned limits, equal, the first fold 
will have a contact with the second fold ; and if C fall below this limit, the 
roots become imaginary, and the two folds become one. Again, if C is 
diminished to the limit where the second equation has the mentioned pair of 
roots equal, the first fold will have a contact with the third ; and when C is 
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less than this, these two folds form but one. And when C is less than both 
these limits, there will be but one fold to the surface. 

In the spaces inclosed by the first and second folds the velocity, relative 
to the moving axes of coordinates, is real ; but, in the space lying between 
these folds and the third fold, it is imaginary ; without the third fold it is 
again real. Thus, in those cases, where C and y have such values that the 
three folds exist, if the body, whose motion is considered, is found at any 
time within the first fold, it must forever remain within it, and its radius 
vector will have a superior limit. If it be found within the second fold, the 
same thing is true, but the radius vector will have an inferior as well as a 
superior limit. And if it be found without the third fold, it must forever 
remain without, and its radius vector will have an inferior limit. 

Applying this theory to our satellite, we see that it is actually within the 
first fold, and consequently must always remain there, and its distance from 
the earth can never exceed 109.694 equatorial radii. Thus, the eccentricity 
of the earth's orbit being neglected, we have a rigorous demonstration of a 
superior limit to the radius vector of the moon. 

In the cases, where C and y have such values that the surface forms but 
one fold, Jacobi's integral does not afford any limits to the radius vector. 

When we neglect the solar parallax and the lunar inclination, the pre- 
ceding investigation is reduced to much simpler terms. The surface then 
degenerates into a plane curve, whose equation, of the sixth degree, is 



li 



— C—ln' 2 x 



It is evidently symmetrical with respect to both axes of coordinates, and is 
contained between the two right lines, whose equations are 



3C — ~f~ -* / . . 



3»' : 

and which are asymptotic to it. It intersects the axis of y, at two points, 
whose coordinates are 



The cubic equation, 



9 = ± G 



r r 4- —-- = 

3n' 2 ^ 3n' 2 
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gives the values of r, for which the curves intersect the axis of x. If 

(2C)§- > 9(in', 
this equation has two real roots between the limits and + \ - —-. If 

(2C)i = %«', 
these roots become equal. And if 

(2tf)S < V, 
there are no real roots between these limits, and the curve has no intersection 
with the axis of x. The figures below exhibit the three varieties of this curve. 



FlS-l 




Fig. Z 



f;&2 




Fig. 1 represents the form of the curve in the case of our moon. In Fig. 2 
we see that the small oval of Fig. 1 has enlarged and elongated itself so as to 
touch the two infinite branches; while, in Fig. 3, it has disappeared, the 
portions of the curve, lying on either side of the axis of x, having lifted 
themselves away from it, and the angles having become rounded off. In Fig. 
1, the velocity is real within the oval, and also without the infinite branches, 
but it is imaginary in the portion of the plane lying between the oval and 
these branches. Hence, if the body be found, at any time, within the oval, 
it cannot escape thence, and its radius vector will have a superior limit ; and, 
if it be found in one of the spaces on the concave side of the infinite branches, 
it cannot remove to the other, and its radius vector will have an inferior 
limit. 
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In the case represented in Fig-. 2, the same things are true, but it seems 
as if the body might escape from the oval to the infinite spaces, or vice versa, 
at the points where the curve intersects the axis of x. However, at these 
points, the force, no less than the velocity, is reduced to zero. For the 
distance of these points from the origin is the positive root of the equation 

or 

V2ff _ r%m' 
3n' ~~ Sn' ' 
and this value is the same as that given by the equation 

JL _ sn' 2 = 0. 

In consequence the forces vanish at these two points, and thus we have two 
particular solutions of our differential equations* 

In the case represented in Fig. 3, the integral does not afford any 
superior or inferior limit to the radius vector. 

The surface, or, in the more simple case, the plane curve, we have 
discussed, is the locus of zero velocity ; and the surface or plane curve, upon 
which the velocity has a definite value, is precisely of the same character and 
has a similar equation. It is only necessary to suppose that the C of the 
preceding formulae is augmented by half the square of the value attributed 
to the velocity. Thus, in the case of our moon, it is plain the curves of equal 
velocity will form a series of ovals surrounding the origin, and approaching 
it, and becoming more nearly circular as the velocity increases. 

Applying the simple formulas, where the solar parallax is neglected, to 
the moon, we find that the distance of the asymptotic lines, from the origin, is 



4. 



= 500.4992. 
Sn' 2 



The distance of the points on the axis of x, at which the moon would remain 
stationary with respect to the sun, is 

ZljL = 235.5971. 

\q„,2 



'Sn'- 



*The corresponding solution, in the more general problem of three bodies, may be seen in the 
Mecanique Celeste, Tom. IV, p. 310. 

7 
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If the auxiliary angle is derived from the equation 

si " • = (H- 

we get 

6 = 32° 49' 6".63 ; 

and the distances from the origin, at which the curve, of zero velocity inter- 
sects the axis of x, are given by the two expressions 

2V2<7 . 

— K : sin — , 

3»' 3 

2V2C . ( aao 6\ 

Sn'- sm V 60 - y) ' 

and the numbers are 109.6772 and 435:5623. These values differ but little 
from the previous more general determinations. 

(To be continued.) 



